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^ . ABSTRACT 

' We report on timing and single-pulse observations of 17 pulsars discovered at 

£S) \ the Arecibo observatory. The highlights of our sample are the recycled pulsars 

Cn ■ J1829+2456, J1944+0907 and the drifting subpulses observed in PSR J0815+0939. 

00 ' For the double neutron star binary J1829+2456, in addition to improving upon our 

existing measurement of relativistic periastron advance, we have now measured the 
*n . pulsar's spin period derivative. This new result sets an upper limit on the transverse 

speed of 120 km s _1 and a lower limit on the characteristic age of 12.4 Gyr. From 
r"| . our measurement of proper motion of the isolated 5.2-ms pulsar J1944+0907, we in- 

fer a transverse speed of 188 ± 65 km s -1 . This is higher than that of any other 
isolated millisecond pulsar. An estimate of the speed, using interstellar scintillation, 
of 235 ± 45 km s _1 indicates that the scattering medium along the line of sight is 
non-uniform. We discuss the drifting subpulses detected from three pulsars in the 
, sample, in particular the remarkable drifting subpulse properties of the 645-ms pul- 

sar J0815+0939. Drifting is observed in all four components of the pulse profile, with 
the sense of drift varying among the different components. This unusual 'bi-drifting" 
behaviour challenges standard explanations of the drifting subpulse phenomenon. 
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■ Key words: pulsars: general — pulsars: individual (PSR J0815+0939, 

PSR J1829+2456, PSR J1944+0907) — radiation mechanisms: non-thermal 



1 INTRODUCTION 

Since the first systematic pulsar survey using the Arecibo 
telescope jHulse fc Tavlorl Il974 Il975allbh which discov- 
ered the original doub le neutron star binary B1913+16 
jHulse fc Tavlorl Il975alh there have been a number of 
searches carried out with this s ensitive system (e.g. 
IStokes et all Il986l : iNice et alJ Il995l) . Timing of the pul- 
sars discovered in these searches has provided a wealth 
of information, including the fir st evidence for the exis- 
tence of gravita tional radiation jTavl or fc Weisberd Il982t 
IWeisberg fc Tavlorlll984D a nd th e first extra-solar planetary 
system dWolszczan fc FraijFl99^ . 

Between 1994 and 1998 the telescope underwent a ma- 
jor upgrade, during which the receivers were parked on 
the meridian but observations could still be made as the 
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sky drifted past. Pulsar surveys conducted in this drift- 
scan mode have resulted in the discovery of 130 pul- 



sars jRav et alJll995l ICamilo et al.l 11996 


: iRav et alJll996l: 


iLommen et alJl2000 


; Lorimer et al. 2004; 


McLaughlin et alJ 


2002:lLorimer et al 


2005). 



In this paper we report on timing and single-pulse ob- 
servations of 17 pulsars discovered at Arecibo. Nine of the 
pulsars were discovered in our recent drift-sc an searches 
jMcLaughlin et aljEool IChampion et alJl2004l) . while the 
remainder were found in earlier surveys and have not 
been subje cted to a follow-up campaign of timing ob- 
servations jRav et alJ Il996t IZeoka et alJll996l: iNicd Il999l: 
iMcLaughlin et al J 12000 ) . We discuss the timing results in- 
cluding the proper motion of millisecond pulsar (MSP) 
J1944+0907 and an updated timing solution for the rela- 
tivistic binary J1829+2456. We then describe the drifting 
subpulses seen in three of the pulsars, including J0815+0939. 
This unique pulsar shows drifting in all four components of 
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its unusual pulse profile, with the subpulses in one compo- 
nent drifting in a different sense to the others. 



2 OBSERVATIONS AND ANALYSIS 

The observations presented here were carried out between 
February f 999 and March 2005 using the Arecibo telescope. 
Most of the observations used the 430-MHz Gregorian dome 
receiver, which has an average gain of 11 K Jy _1 and an av- 
erage system temperature of 45 K, in combination with the 
Penn State Pulsar Machine (PSPM) , a 1 28-channel analogue 
filterbank spectrometer dCadwelll 119971) . During the early 
stages of timing when the periods of the pulsars were not 
well determined, the PSPM was used in search mode. The 
PSPM samples incoming voltages from the telescope every 
80/xs over a bandwidth of 7.68 MHz. These data were then 
dedispersed by delaying successive channels in time corre- 
sponding to the nominal dispersion measure (DM) of each 
pulsar to produce time series which were subsequently folded 
to produce integrated pulse profiles. To confirm that the pe- 
riod was the fundamental rather than a harmonic, the data 
was folded at twice the period and the profile was checked. 
Once the period of each pulsar had been determined ac- 
curately, the timing mode of the PSPM was used. In tim- 
ing mode the data in each frequency channel are folded on- 
line. The 128 individual profiles can then be dedispersed and 
summed to produce a single profile. 

The first stage of timing analysis is to produce a high 
signal-to-noise ratio (S/N) profile to use as a template. Ini- 
tially this was done by shifting the profile so the peak was 
at a phase of 0.5 before summing the profiles. The time of 
arrival (TOA) of each profile was then calculated by con- 
volving the profile with the template and adding this offset 
to the start time of the observation. These TOAs were then 
analysed with the TEMPO software package to produce a 
preliminary phase-connected solution (where every rotation 
of the pulsar is accounted for) over the time-span of the ob- 
servations. In most cases a simple spin down model with five 
free parameters (pulse phase, period, period derivative, right 
ascension, and declination) was used. For PSR J1944+0907 
the proper motion is a significant effect and is included in 
the fit. The relativistic binary pulsar J1829+2456 requires 
an additional five classical and two relativistic parameters 
to produce a phase connected solution. For each pulsar the 
resulting ephemeris was used to predict the pulse phase 
of each profile. These were then aligned before they were 
summed to produce a better template which was subse- 
quently used to generate new TOAs, weighted by S/N, which 
were analysed using TEMPO. This additional iteration of- 
ten resulted in significant reduction of the root-mean-square 
observed minus model timing residual. With the exception 
of PSR J1917+0834, a single template was sufficient for the 
timing analysis. For J1917+0834, due to the short length 
of the integrations and long period of the pulsar, only a 
small number of pulses were summed causing the relative 
strengths of the two peaks in the profile to vary from obser- 
vation to observation, a template for each mode was required 
to provide accurate TOAs in the initial stages. 

1 http:/ /pulsar. princeton.edu/tempo 



Most of the previously published DM measurements 
were based on the original observations and so were poorly 
constrained. In December 2004 the pulsars were observed at 
320, 327, 334 and 430 MHz using the PSPM. The TOAs 
produced were analysed using TEMPO by holding all the, 
now well constrained, parameters of the ephemeris constant 
and fitting for DM. This new value of DM was then used to 
dedisperse the data again and improve the template, TOAs 
and ultimately the ephemerides. 

To calibrate the intensity scale of the pulse profiles, 
we followed the procedure described by Lorimer, Camilo 
& Xilouris (2002) and scaled each profile by the factor 
1000T/ (DG) mjy, where T is the total system noise tem- 
perature (K), D is the DC offset of the profile (in PSPM 
counts) and G is the antenna gain (K Jy _1 ). Due to the 
design of the Arecibo telescope, both T and G are strong 
functions of zenith angle. This dependence is described an- 
alytically by Perillat (1999). The total system temperature 
also includes the bac kground sky temperature contribution 
jHaslam et alj Fl9821. The error is based on the standard 
deviation of the observed fluxes. 



3 TIMING RESULTS 

The basic timing parameters, pulse widths and flux densi- 
ties are summarised in Table Q Listed are the pulsar name, 
position (in right ascension and declination), period P, its 
derivative P and epoch, DM, span of observations, number 
of TOAs and the post-fit weighted rms of the timing residu- 
als. Table|5|lists the additional timing parameters for the bi- 
nary pulsar J 1829+2456 and the isolated MSP J 1944+0907. 
The 430-MHz profiles derived from a phase-coherent sum 
over many individual observations are shown in Fig. 

Further observed and derived parameters are sum- 
marised in Table |3] Listed are the pulse widths at 10%, 
Wio, and 50%, logo, of peak amplitude; the equivalent pulse 
width (the width of a top hat function with the same inte- 
grated flux and peak amplitude), w eq ; the mean 430-MHz 
phase-averaged flux density, S430; distance d (inferred from 
the free electron density model of Cordes & Lazio 2002); 
Galactic height [z = dsinfe, where b is the Galactic lati- 
tude); 430-MHz luminosity (L430 = 5*43od 2 ); spin-down age 
(r = P/2P); spin-down energy-loss rate (E = An 2 IP/P 3 , 
where the neutron star moment of inertia / is taken to 
be 10 45 g cm 2 ); and inferred surface dipole magnetic field 

strength (B = 3.2 x W 19 \/p~P G; Lorimer & Kramer 2005). 

In Fig. |5|we show the locations of the 17 pulsars on the 
period-period derivative (P-P) diagram. As can be seen, the 
majority of the pulsars occupy the lower end of the 'island' 
of points suggesting that they are relatively old members 
of the normal population. This is in keeping with previ- 
ous samples of pulsars d rawn from low-frequency surveys 
(e.g. lLorimer et all 12002) which are less effective at find- 
ing young pulsars than surveys at higher frequencies (e.g. 
iManchester et al-lEoOlh . 

3.1 PSR J1829+2456 

This 41-ms pulsar is in a mildly eccentric binary sys- 
tem where the compa nion is most likely a neutron star 
JChampion et alJl2004h . The rate of advance of periastron, 
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Table 1. Observed parameters 



Name R.A. (J2000) Docl. (J2000) P P Epoch DM Span IV TOAs RMS 

PSR (h:m:s) (":':") (s) (10~ 15 ) (MJD) (pc cm -3 ) (MJD) (is 



J0152+0948 


01:52:23.73(5) 


+09:48:10(2) 


2.74664729014(10) 


1.700(4) 


52858 


21.87(2) 


52597- 


-53004 


26 


350 


J0546+2441 


05:46:28.76(3) 


+24:41:21(6) 


2.84385038524(6) 


7.65(3) 


52914 


73.81(3) 


52707- 


-53119 


29 


706 


J0815+0939 


08:15:08.776(9) 


+09:39:50.7(6) 


0.64516119106(5) 


0.139(6) 


52854 


52.667(9) 


52600- 


-52796 


53 


270 


J1503+2111 


15:03:54.60(4) 


+ 21:11:09.3(5) 


3.31400150122(12) 


0.14(2) 


52857 


11.75(6) 


52595- 


-53118 


12 


941 


J1746+2245 


17:46:00.84(3) 


+22:45:28.9(8) 


3.46503778264(6) 


4.919(7) 


52765 


52.0(20) 


52401- 


-53129 


19 


1240 


J1807+0756 


18:07:51.235(2) 


+07:56:43.39(5) 


0.464300493199(2) 


0.1294(2) 


53064 


89.29(3) 


52765- 


-53362 


22 


134 


J1821+1715 


18:21:13.503(4) 


+ 17:15:47.06(7) 


1.366682059422(10) 


0.8714(11) 


53064 


60.47(7) 


52765- 


-53362 


23 


240 


J1829+2456 


18:29:34.6668(2) 


+ 24:56:18.193(3) 


0.04100982358117(4) 


0.0000525(15) 


52887 


13.8(5) 


52785- 


-53447 


50 


22 


J1843+2024 


18:43:26.18(7) 


+20:24:54.6(12) 


3.4065386707(4) 


1.04(4) 


53064 


85.3(2) 


52765- 


-53362 


19 


3829 


J1849+2423 


18:49:34.6763(4 ) 


+24:23:45.976(12) 


0.27564149867101(7) 


0.152582(4) 


52199 


62.239(5) 


51194- 


-52923 


244 


409 


J1915+0752 


19:15:01.94(2) 


+07:52:09.2(6) 


2.05831378861(7) 


0.139(9) 


53048 


105.3(3) 


52733- 


-53362 


60 


2891 


J1916+0748 


19:16:51.5(2) 


+07:48:00(5) 


0.54175211381(17) 


10.708(17) 


53048 


304.0(15) 


52733- 


-53362 


15 


14754 


J1917+0834 


19:17:48.853(6) 


+08:34:54.63(14) 


2.129665364024(4) 


17.4935(2) 


52345 


29.18(6) 


51327- 


-53362 


36 


608 


J1933+0758 


19:33:19.407(2) 


+07:58:07.01(7) 


0.437454407217(2) 


0.2183(2) 


53048 


165.02(4) 


52732- 


-53362 


19 


171 


J1944+0907 


19:44:09.32084(3) 


+09:07:23.2388(12) 


0.0051852019041260(3) 


0.00001713(2) 


52913 


24.34(2) 


52732- 


-53447 


27 


2.6 


J2007+0910 


20:07:58.1031(12) 


+09:10:12.93(6) 


0.4587348258074(6) 


0.33244(12) 


52994 


48.6(2) 


52595- 


-53392 


25 


116 


J2045+0912 


20:45:47.3438(9) 


+09:12:29.16(4) 


0.3955551117394(5) 


0.19520(7) 


52997 


31.776(4) 


52601- 


-53392 


25 


78 



Observed parameters of the 17 pulsars derived from the TEMPO timing analysis described in section 2. Figures in parentheses are 1-cr 
uncertainties in the least significant digits. The arrival times from which these ephemerides are derived are freely available online as 
part of the European Pulsar Network (EPN) database jhttp: //www .jb.man.ac.uk/~pulsar/Resources/epnf . 
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Figure 1. Integrated 430-MHz pulse profiles. All profiles show full pulse phase over 256 bins with the exception of PSR J1944+0907 
which has 128 bins. PSR J1944+0907 is dispersion-smeared over ~4 bins. These profiles are freely available as part of the EPN database. 



uj, is now measured to a significance of 180-er. Interpreting uj 
as purely relativistic precession implies that the total mass 
of the binary M = 2.59 ± 0.02 M (see Fig. BJ. The lower 
companion mass limit of 1.26 M , and upper limit to the 
pulsar mass of 1.34 Mq are given by the constraint that the 
inclination angle i ^ 90° . A very conservative upper limit on 
the companion mass M c = 2.61 Mq , and a minimum inclina- 
tion angle i > 29° , can also be inferred. A more likely upper 
limit of M c = 1.33 Mq is obtained by assuming a minimum 
pulsar mass of 1.26 Mq (the lowest mea sured mass of an y 
neutron star, namely PSR J0737-3039B: iLvne et ai]|2004) . 
An upper limit of 8 ms can now also be placed on the rel- 
ativistic time-dilation and gravitational redshift parameter, 
7, but, as shown in Fig. 01 this is not yet constraining. The 
determination of these limits is described in more detail in 
IChampion et alj (120041) . 



The extended baseline of our TOAs has also resulted in 
a measurement of P for this pulsar. However this P should 
be treated with caution since any transverse component of 
proper motion of the system would lead to secular accel- 
eration dShklovskiill970t ICamilo et ailll994l) and will con- 
tribute to the observed period derivative, 

i = lY f = ^*l^(f)(^^)\ (i) 
P c a \ kpey \masyr -1 J 

where P contribution due to the proper motion, Vr is the 
transverse speed and /xt is the proper motion. Indeed, as- 
suming the current measurement of P is entirely due to 
proper motion, and adopting the nominal DM-derived dis- 
tance of 1.2 kpc, the Shklovskii effect places an upper limit 
of 20 mas yr _1 on the proper motion, corresponding to a 
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Table 2. Additional parameters for pulsars J1829+2456 and J1944+0907 



PSR 



J1829+2456 




x, Projected semi-major axis (lt-sec) 


7.238(2) 


Pi,, Binary period (days) 


1.176027941(16) 


e, Eccentricity 


0.1391412(19) 


uj, Longitude of periastron (deg) 


229.92(2) 


To, Epoch of periastron (MJD) 


52848.579775(3) 


u), Rate of advance of periastron (deg yr — 1 ) 


0.2919(16) 


7, Relativistic time-dilation and gravitational redshift (ms) 


3.6(38) 


Mean TOA uncertainty (fis) 


33.7 


J1944+0907 




Ha, Proper motion in R.A. (mas yr —1 ) 


12.0(7) 


Proper motion in Decl. (mas yr — 1 ) 


-18(3) 


Mean TOA uncertainty (/is) 


3.6 



Figures in parentheses are 1-cr uncertainties in the least significant digits. 



Table 3. Further observed and derived parameters 



Name w w w 50 w eq S430 d z L430 Log t Log E Log B 

PSR (ms) (ms) (ms) (mjy) (kpc) (kpc) (mjy kpc 2 ) (kyr) (ergs s — 1 ) (G) 



J0152+0948 


103 


11 


17 


0.91(3) 


0.9 


-0.71 


0.7 


4.4 


30.5 


12.3 


J0546+2441 


51 


28 


25 


2.65(10) 


2.0 


-0.07 


10.6 


3.8 


31.1 


12.7 


J0815+0939 


183 


166 


75 


3.7(2) 


2.5 


0.98 


23.1 


4.9 


31.3 


11.5 


J1503+2111 


80 


26 


39 


1.3(2) 


1.0 


0.88 


1.3 


5.6 


29.2 


11.8 


J1746+2245 


125 


34 


50 


1.26(7) 


3.0 


1.22 


10.8 


4.1 


30.7 


12.6 


J1807+0756 


27.3 


18.2 


17.6 


2.16(6) 


3.5 


0.81 


26.5 


4.8 


31.7 


11.4 


J1821+1715 


64 


45 


39 


3.7(2) 


2.8 


0.70 


29.0 


4.1 


31.1 


12.0 


J1829+2456 


1.7 


0.8 


1.0 


0.3(1) 


1.2 


0.33 


0.4 


7.1 


31.5 


9.2 


J1843+2024 


29 


40 


29 


0.37(3) 


4.0 


0.75 


5.9 


4.7 


30.0 


12.3 


J1849+2423 


26.1 


14.4 


14.6 


2.3(6) 


3.2 


0.63 


23.6 


4.5 


32.5 


11.3 


J1915+0752 


60 


24 


28 


1.70(10) 


3.8 


-0.11 


24.6 


5.4 


29.8 


11.7 


J1916+0748 


230 


80 


50 


1.14(9) 


8.1 


-0.29 


74.8 


2.9 


33.4 


12.4 


J1917+0834 


81 


68 


35 


0.44(10) 


1.9 


-0.06 


1.6 


3.3 


31.9 


12.8 


J1933+0758 


14.6 


6.8 


6.8 


0.34(3) 


6.0 


-0.58 


12.2 


4.5 


32.0 


11.5 


J1944+0907 


2.5 


0.5 


0.5 


3.9(3) 


1.8 


-0.23 


12.6 


7.1 


33.2 


8.3 


J2007+0910 


16.6 


8.3 


8.7 


1.51(5) 


2.0 


-0.55 


6.0 


4.3 


32.1 


11.6 


J2045+0912 


25.5 


10.4 


13.7 


3.5(2) 


1.9 


-0.66 


12.6 


4.5 


32.1 


11.5 



Pulse profile widths, 430-MHz flux density measurements and various derived parameters for the 17 pulsars. All distance measurements 

have statistical error of 25% as estimated by Cordes & Lazio (2002), individual errors may be larger. Figures in parentheses are 1-cr 
uncertainties in the least significant digits for the observed parameters. For PSR J1944+0907, the Shklovskii-corrected P has been used 
to calculate the luminosity, spin-down age, spin-down energy loss rate and inferred surface dipole magnetic field strength. 

transverse speed of ~118 km s -1 . This transverse speed is 
consistent with other double neutron star systems so it is 
possible that the intrinsic P is considerably lower. 

However if the P of PSR J1829+2456 is assumed to be 
intrinsic it places this double neutron star system in an area 
of the P-P diagram inhabited by various eccentric binary 
systems (see Fig- EJ. This P implies a characteristic age of 
at least 12.4 Gyr and a magnetic field strength of 1.5 x 10 9 
G. As for other recycled pulsars this large characteristic age 
suggests that its birth spin period was comparable to its 
current spin period. Of the other doub le neutron star sys - 
tems, only J1518+4904 has a lower P jHobbs et alJl2004l) . 
PSR J1829+2456 is currently the subject of an extended 
timing campaign to measure a proper motion and so remove 
any contribution of the Shklovskii effect from the measure- 
ment of P. Simulated TOAs suggest that a proper motion 
of 15 mas yr _1 should be measured to ~ 30% and the ex- 



pected 7 (1.3 ms) should be measured to ~ 40% accuracy 
within the next year. Simulations also suggest that it would 
be possible to measure the Shapiro delay with the current 
observing equipment if the inclination was ~ 90°. 

3.2 PSR J1944+0907 

PSR J1 944+0907 is an iso l ated 5. 2-ms pulsar with a broad 
profile jMcLauehlin et alj I2005I) . Isolated MSPs are rare 
objects; only 13 are currently known jLorimerll2005f) . While 
binary MSPs are thought to be spun up by the accretion of 
matter from an evolved binary companion, the formation 
process for isolated MSPs is presently not clear. Either 
they form through a different formation channel (such as 
accretion- induced collapse), or they ablate their companio n 
with their strong relativistic wind iRuderman et alJ ll989L 
As shown in Fig. 111 the proper motion of this pulsar, llt, is 
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Figure 2. The P — P diagram showing all currently 
known pulsars listed in the on-line Australia Tele- 
scope National Facility (ATNF) pulsar calatalogue 
(http://www.atnf.csiro.au/rcsearch/pulsar/psrcatl. The stars 
indicate the 17 pulsars in the current sample. Young pulsars 
are in the top right of the plot while old recycled pulsars are 
in the bottom left. The line indicates the age of th e Universe, 
r = 14 Gyr, equivalent to Hq = 70 km s _1 Mpc -1 J York et alj 
120051) . The P of PSR J1944+0907 has been Shklovskii-corrected 
assuming the DM-inferred distance (see section 13.21 . The P of 
PSR J1829+2456, for which there currently no proper motion 
measurement, has not been Shklovskii-corrected (see section 

rm . 



detected very significantly in our timing observations from 
which we obtain fii = 22 ± 2 mas yr _1 . This would imply 
a transverse speed Vr = /J-rd = 188 ± 65 km s _1 ; this is 
significantly higher than the average for recycled systems 
(87 ± 13, Hobbs et al. 2005) and, with the exception of 
PSR B1257+12 which has a planetary system, this is by 



far the highest of all isolated 


and indeed most binary) 


MSPs so far JCamilo et al.lll993L 


1996 


Nice & Taylor 1995; 


Coenard et alJll995l Y] 


?oscano et alJ|l999: Wolszczan et alJ 


2000l: lLommen et al. 


l2000t 


Edwards & Bailcs 2001; 


Lorimer et alJ 12004: 


Faulkner et alJ 


2004; Hobbs ct al. 


2004). 



The very large transverse speed implied by the proper 
motion would have a significant effect on the observed P as 
described in 13.11 When this Shklovskii effect is taken into 
account the intrinsic P drops to 6.2 x 10 -21 . This value 
has been used to calculate the derived parameters listed in 
Table 01 

PSR J1944+0907 also scintillates strongly at 430 MHz, 
as shown in Fig. An auto-correlation analysis of its dy- 
namic spectrum indicates a bandwidth and timescale of Avd 
= 90 kHz and Atd = 101 s. The interstellar scintillation de- 
rived transverse speed, Viss, is given by 



Viss = 



AdiAvJ 



vAt d 



(2) 



where the constant A — 2.53 x 10 4 km s _1 for a uni- 
form Kolmogorov scattering medium, d is the distance in 
kpc, v is frequency in GHz, Avd is scintillation band- 
width in MHz and Atd is the scintillation timescale in 
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Figure 3. Current constraints on the masses of PSR J1829+2456 
and companion. The grey bar shows the mass limits for neu- 
tron stars, the upper limit described by Thorsett & Chakrabarty 
(1999), and the low measured mass of PSR J0737-3039B (Lyne 
et al. 2004). There is currently only an upper limit for the time- 
dilation and redshift parameter, 7. The upper and lower limits of 
ill are now so close they appear as one thick line. The masses of 
the system are constrained to be between the two w limits and 
above the i = 90° limit. The marginal figures show the probability 
density functions for the pulsar and companion masses assuming 
a random distribution of inclination angles. The lines of constant 
mass indicate the median of the distributions. 
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Figure 4. The timing residuals for PSR J1944+0907 before fit- 
ting for proper motion. The solid line indicates the fit for a proper 
motion of 22 mas yr — 1 . 



seconds dCordes fc Rickettlll998l) . Assuming a distance of 
1.8±0.5 kpc, the implied transverse speed is 235T45 km s _1 . 

Because the transverse speeds inferred from proper mo- 
tion and scintillation depend differently on the distance, we 
may use these measurements to calculate a distance consis- 
tent with both measurements. Assuming a uniform scatter- 
ing medium, this distance is 3.6 kpc, with implied transverse 
speed of 340 km s _1 . Not only is this far higher than the 
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428 430 432 
Frequency (MHz) 

Figure 5. Dynamic spectrum (i.e. on-pulse intensity versus time 
and frequency) for J1944+0907 at 430 MHz. The intensity is in- 
dicated by the linear greyscale (black: low intensity, white: high 
intensity). The spectrum was formed with 10-s integrations us- 
ing the PSPM, with 128 frequency channels across the 7.68-MHz 
bandwidth. 

transverse speeds measured for any isolated MSPs but would 
imply a Shklovskii correction greater than the observed P! 
We therefore conclude that the scattering medium is not 
uniform and that the actual distance and speed are closer 
to 1.8 kpc and 188 km s _1 , respectively. 



4 DRIFTING SUBPULSE BEHAVIOUR 

Single pulses often consist of several components, called 
subpulses, which change position, often drifting from pulse 
to pulse. This phenomenon was first reported by Drake 
& Craft (1968) and has since been recognised in roughly 
100 radio pulsars (Weltevrede, Edwards & Stappers 2005). 
It is usually understood w i thin the E X B drift model 
of iRuderman fc Sutherland! jl975) in which radio emission 
comes from discrete locations, or "sub-beams", positioned 
on a circle around the magnetic pole. The sub-beams rotate 
and move in and out of our line of sight. The Ruderman 
& Sutherland model has been fairly successful in explain- 
ing pulsar drifting patterns . In fact, for PSR B0943+10, 
iDeshpande fc Rankinl jl999T) have been able to determine 
that there are 20 sub-beams rotating around this pulsar's 
magnetic axis, with a rotation timescale close to 100 s. 

In some cases our PSPM search-mode data were of suf- 
ficient quality to analyse individual pulses. In our sample 
three of the pulsars show drifting subpulses. The drifting 
can be c haracterised by the drift rate and two periods (P2 
and P3; lBackerlll973l) . The drift rate is simply the longi- 
tude of pulse phase drifted by each subpulse per second. 
The separation in longitude between consecutive subpulses 
and the separation at a given longitude between successive 
subpulses are denoted by P2 and P3. The single pulses for 
the three pulsars are shown in Fig. ©and the measured drift 
rate and P3 for each pulsar are listed in Table [I] P2 is not 
applicable in these cases as there is never more than one 
subpulse (per component) visible during a single pulse. For 
PSR J0815+0939, P3 was determined using a fluctuation 



spectrum analysis, in which a Fourier transform is applied 
to a series of single pulses (in this case 256) for a range of 
pulse longitudes. 

The single-pulse behaviour of the 645-ms pulsar 
J0815+0939 is remarkable. This pulsar's profile is unusual, 
consisting of four components with separations of 41, 58, and 
54 ms and not easily desc ribed by stand ard phenomenolog- 
ical classification schemes iRankinl ill 9831) . As shown in Fig. 
HJ drifting is observed within all four components. We refer 
to these (from left to right) as components I, II, III and IV. 
The drift within component I is difficult to measure and ap- 
pears to change sense within an observation. The drift rates 
of components II, III and IV are very stable from epoch to 
epoch and do not appear to change sense. However, the drift 
sense of component II differs from those of components III 
and IV. While of opposite sign, the drift rates of components 
II and IV are equal within the errors, indicating a possible 
relationship. Measurements at 327 and 430 MHz indicate an 
evolution of the drift rate with frequency, with a higher drift 
rate at the lower frequency. 

While drift sense reversals within a single component 
are quite common and are attributed to aliasing effects 
jGupta et alJl2004l) . it is rare for different components to 
show different drift senses (Hankins & Wolszczan 1987; 
Weltevrede et al. 2005). How do we interpret this strange 
behaviour? It is very difficult to reconcile with the stan- 
dard Ruderman & Sutherland model in which sparks travel 
around a circular region centred on the pulsar's magnetic 
axis. In this case the four components would probably be 
associated with two nested hollow-cone beams. Similar drift 
sense would be expected in components I and IV and also II 
and III but this is not observed. This pulsar's " bi-dri fting" 
behav iour has been geometrically modelled by iQiao et alJ 
i2004l) . While promising, this model is based on one short 
integration at one frequency. Edwards & Stappers (2003) in- 
terpret the behaviour of B1918+19 as being due to multiple 
imaging of subbeam systems. Further studies are planned 
which will determine if there are any aliasing effects while 
polarisation measurements will help determine the relation- 
ship between the multiple components. 

In addition to J0815+0939, both PSR J1916+0748 and 
PSR J2007+0912 show weak drifting as shown in Fig.© For 
J1916+0748, the broad integrated profile can be explained 
as the superposition of a series of drifting subpulses. Further 
studies of these pulsars with higher sensitivity and time res- 
olution in future may well be worthwhile. 



5 CONCLUSIONS 

We have presented a timing and single-pulse study of a 
sample of 17 pulsars discovered with the Arecibo telescope. 
These observations provide new insights into a wide range 
of neutron star physics. We now briefly summarize the new 
results obtained from the present study. 

Extending the timing baseline of the relativistic binary 
pulsar J1829+2456 presented by Champion et al. (2004) 
has allowed the mass of the system to be further con- 
strained, M = 2.59 ± 0.02 Mq, and an upper limit to be 
placed on the relativistic time-dilation and gravitational 
redshift parameter, 7 < 8 ms. The measurement of P = 
5.25 ±0.15 x 10~ 20 has placed an upper limit on the proper 
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Table 4. Parameters obtained from our analyses of the single-pulse observations described in Section FT! 
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Name 


Frequency (MHz) 


Component 


Drift rate (deg s 1 ) 


ftOO 


J0815+0939 


430 


I 




10.87(16)" 






11 


-1.60(13) 


10.87(16)° 






111 


2.46(19) 


10.87(16)" 






IV 


1.57(16) 


10.87(16)" 




327 


I 




8.72(16)° 






11 


-2.3(5) 


8.72(16)° 






111 


3.2(8) 


8.72(16)° 






IV 


2.68(18) 


8.72(16)° 


J1916+0748 


430 


I 


5.4(14) 


8.6(8) 


J2007+0912 


430 


1 


1.7(6) 


5.4(7) 


Figures 


in parentheses are la- 


uncertainties 


in the least significant digits. 



These parameters are derived from the highest S/N observation for each pulsar, but the drifting properties appear to be consistent 
from epoch to epoch. A more detailed study of the time evolution of the drifting properties of J0815+ 0939 will be presented elsewhere, 
a: Determined from fluctuation spectra. Uncertainties are determined by the resolution of the spectra. 



J0815+0939 (430 MHz) J0815+0939 (327 MHz) J1916+0748 (430 MHz) J2007+0910 (430 MHz) 




Pulse Phase (degrees) Pulse Phase (degrees) Pulse Phase (degrees) Pulse Phase (degrees) 

Figure 6. The single pulses for pulsars J0815+0939, J1916+0748 and J2007+0912. The intensity is indicated by the linear greyscale 
(white: low intensity, black: high intensity). The drifting can be seen as diagonal stripes. In PSR J0815+0939 it is clear that the second 
component drifts in the opposite sense to the third and fourth. Drifting in the first component is unclear and seems to show some sense 
reversal. 



motion, /j,t =20 mas yr -1 , (using the Shklovskii effect) and 
a lower limit on the characteristic age, r ~ 12.4 Gyr. Further 
timing will allow a measurement of this pulsar's proper mo- 
tion and a determination of the intrinsic P. A measurement 
of the parallax through timing is less likely as the current 
level of precision. 



The proper motion of PSR J1944+0939 has been mea- 
sured and revealed the highest transverse speed of any iso- 
lated MSP, V T = 188 ± 65 km s" 1 . An alternative esti- 
mate of the speed using interstellar scintillation, Viss = 
235 ± 45 km s _1 , demonstrates that the scattering medium 
along the line of sight is non-uniform. These measurements 
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suggest that the space velocity of J1944+0939 is higher than 
for any other isolated MSP. Further timing observations will 
ultimately allow the parallax of this pulsar to be measured 
which will resolve current uncertainties in the transverse 
speed estimates. 

The drifting observed in PSR 0815+0939 is very un- 
usual and warrants further study. All four components of 
the pulse profile show drifting subpulse behaviour. Of the 
three components in which we see no drift sense reversals, 
one of the components has a different drift sense to the oth- 
ers. This is difficult to explain using the standard Ruder- 
man & Sutherland model. Further observations, including 
polarimetry, are required to shed light on this remarkable 
pulsar and the implications it may have for emission mech- 
anism models. 
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